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, can be used in purpose of analysis of other asymptotic series in turbulence. 
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1 Introduction 

Theoretical investigations of two dimensional (2D) cross-field transports in the operating 
ITER-FEAT calls for an increasing confidence in the modelling efforts that force one to search 
for the new principles of simulations. 

It has been reported [H |2] that the solutions of forced dynamical equations describing the 
fluid flows (such as the 3D Fourier transformed Navier-Stokes (ENS) equation, the Burgers 
equation, etc.) may be represented, point wise, as the expected values of a branching random 
walk process statistic. The statistic samples the Eourier transforms of both the initial data 
and the forcing at various random frequencies as directed by the leaves of the binary branching 
process, and then combines the results together, in a multiplicative way, according to the binary 
bifurcation nodes of the process. The branching random walk holds in state > 1, A; G N for 
an exponential time with parameter uk"^ ( z/ > is the diffusion coefficient) then it is replaced 
by no offspring with probability 1/2, otherwise replaced by two frequencies j and k — j, with 
j picked uniformly over [1, . . . k — 1]. The process is repeated independently for each of the 
offspring sequences. 

The use of finite binary tree algorithms is very plausible since they can be easily implemented 
and optimized that would drastically improve the effectiveness of large scale simulations on fluid 
flows. However, it has been discovered that the deflnition of the stochastic process relevant 
to the modelling of complicated cross-fleld transports in plasma is a tough problem. Even 
if such the branching multiplicative process exists, its fast convergence is questionable. To 
illustrate the problem, we have performed a preliminary simulation on the ENS. We have tried 
to reproduce the well-known Kolmogorov's inertial range for the fully developed turbulence in 
the 3D ENS forced in the small moments (in the large scales) . We have sampled the binary trees 
independently for each of velocity components, but limit the total number of samples just to 250 
for each point {k, t). For the small moments, just a few different binary trees can be constructed 
and all of them arise in samplings that flts well the fully developed turbulence statistics (see 
Eig. 1). However the number of different branching trees that could appear in simulations 
is huge for k large and apparently grows very fast with k. In this case, the set of randomly 
generated 250 samples obviously does not provide a reliable base for the turbulent statistics 
that is evident from the excessively strong trembling of data worsening as the momentum k 
grows up. The data is presented in Eig. 1 in the form of a box plot indicating the scattering. 
The central point of a box shows the median, the lower line of a box indicates the flrst quartile, 
and the upper line shows the third quartile. The whiskers are extended to the maximum and 
minimum points. Apart from the forced region, an inertial range appears for any nontrivial 
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initial condition. The strong data scattering in the region of large moments indicates that the 
number of samples should be essentially increased. It is important to mention that the classical 
methods developed in [T] do not provide us with tools for the control of convergence for the 
stochastic multiplicative branching processes that forces us seriously concerns to the problem. 

In the present paper, we propose an alternative method that would optimize the simulations 
on fluid flows. It is closely related to that of branching Brownian motion and uses the techniques 
borrowed from the quantum fleld theory to control the computation accuracy. Applications of 
fleld theory methods to the various problems of statistical physics and critical phenomena 
theory have a long history. Nowadays, they have converted into a powerful tool for the effective 
studying of stochastic phenomena described by the partial differential and integro-differential 
equations. The problems of statistical physics, in their turn, had also enriched the fleld theory 
methods as well. 

We study the Green functions of equations describing the incompressible fluid flows by 
constructing their functional representations. Our approach is based on the ideas of the MSR 
formalism (after Martin-Siggia-Rose) treating the classical Langevin equations as a kind of 
Brownian motion. In such a framework, the stochastic averages with respect to the Gaussian 
distributed random forces stirring the system can be interpreted as some functional averages 
(see Sec. 2). It is worth mentioned that they are nontrivial even if the stochastic force is put 
to zero. This fact allows modifying the initial MSR approach in a way to apply it to the forced 
dynamical equations in which the external forcing is not random. The Green functions can be 
computed by means of perturbation theory in the form of Feynman diagram series (Sec. 3). 

We have demonstrated that each Feynman graph in the series equals to an average over a 
forest of multiplicative branching binary trees (of the certain topological structure) implemented 
in p[] that establishes the direct relation between their methods and the fleld theory approach 
that we use (Sec. 3). However, the diagram series converges much faster than the multiplicative 
branching process since each diagram represents an average of branching process over all allowed 
conflgurations of moments. We have demonstrated the advantage of method for the stirred 
Burgers equation supplied with a periodic boundary condition (Sec. 4). In this simple example, 
the diagrams can be computed analytically, and the series converges very fast. The Fourier 
spectrum for the Green function forms the asymptotic series starting from the standard diffusion 
kernel. Fluctuations arise due to the consequent bifurcations of media resulting in the cascade 
of consequent partitions of moments, k = q+ (k — q). The magnitude of relevant corrections to 
the diffusion spectrum tends to zero as t to, and the saddle-points (instanton) analysis can 
be applied to study the "large order" asymptotic contributions (see Sec. 5). The central point 
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of method is that a coefficient in the asymptotic series is calculated by the Cauchy integral 



where the integration contour encircles g = in the complex plane. If the asymptotic series 
G^^^g^ is a functional integral, but is a large number, one can estimate the large order 
coefficient G^^^ calculating the functional integral by the steepest descent method jH]. The 
saddle-point equations describing the "equilibrium" state of the hydrodynamical system with 
respect to the large - order contributions into the Green function can be solved as t ^ to when 
the bifurcations rise the corrections to the pure relaxation dynamics. 

The calculations show that the asymptotic coefficients demonstrate the factorial growth like 
the most of models in quantum field theory. The asymptotic series for the Green function can 
be summarized by means of the Borel procedure (see Sec. 5). In the limit t —>■ to the corrections 
to the diffusion kernel have the closed analytical form. 

The proposed approach can be broadly used in the many problems of forced dynamics for 
the optimization of numerical computations and the analysis of asymptotic series. We conclude 
in the last section. 

2 Stochastic dynamics as the Brownian motion 

It is well known that many problems in stochastic dynamics can be treated as a generalized 
Brownian motion {6 {u ~ u {x,t)))^, in which the classical random field indicating the position 
of a particle u{x,t) meets a Langevin equation 



where ^ is the Gaussian distributed stochastic force characterized by the correlation function 
= Here the angular brackets (. . .)^ denote an average position of particle with respect 

to the statistics of C,- Q{u) is some nonlinear term which depends on the position u{x,t) and 

its spatial derivatives. Such a representation was a key idea of the formalism 

An elegant way to obtain the field theory representation of stochastic dynamics is given by 

the functional integral 



where the Tr-operation means the integration J dxdt and the summation over the discrete 
indices. The instantaneous positions u {x, t) meet the dynamical equation that can be taken 



2m 



/ ^ E «"G<"' = ipj'it J"G<") exp (-N logs) 



u (x,t) = Q{u) + ^, 




(1) 
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5{u-u{x,t)) ^ S{u{x,t) -Q{u) (2) 

should the solution of dynamic equation exists and is unique. The use of integral representation 
for the 5— function in (Q) transforms it into 



j D^DuDu'exp Tr (^~^D^^ - iiu' + u'Q {u) + u'^ det M, 



(3) 



in which u' (x, t) is the auxiliary field that is not inherent to the original model, but appears 
since we treat its dynamics as a Brownian motion. The Jacobian det M relevant to the change 
of variables Q is discussed later. 

The Gaussian functional integral with respect to the stochastic force ^ in Q is calculated 



j DuDu'exp (S {u, r])) det M, S{u, u') = Tr 



-u'D^u' — u'u + u'Q (u) 



(4) 



By means of that all configurations of ^ compatible with the statistics are taken into account. 
The integral (jH)) identifies the statistical averages (. . .)^ with the functional averages of weight 
exp S. The formal convergence requires the field u to be real and the field u' to be purely 
imaginary. 

The functional averages in (0} can be represented by the standard Feynman diagram series 
exactly matching (diagram by diagram) the usual diagram series found by the direct iterations 
of the Langevin equation averaged with respect to the random force This fact justifies the use 
of functional integrals in stochastic dynamics at least as a convenient language for the proper 
diagram expansions. 

The Jacobian det M in depends upon the nonlinearity Q (u) . If Q (u) does not depend 
upon the time derivatives, all diagrams for det M are the cycles of retarded lines A oc 6 (t — t') 
and equal to zero excepting for the very first term, 

det M = const ■ exp TrA, (5) 

Following in the paper, we take the convention for the Heaviside function of zero argument, 
6(0) = 0, so that detM = const. 

Below, we consider the field theory formulations of equations describing the fluid flows. 
We suppose that the fluid is incompressible, div u = 0. In the Fourier space, the relevant 
transversal (rotational) component of u can be allocated by means of the transversal projector. 
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The transversal projection of Navier-Stocks equation (NS) lacks of the pressure gradient and 
the longitudinal component of nonlinear term, [{ud) Any solution of it corresponds to some 
vortex of momentum k. The MSR action functional for the transversal NS is 



S (u, u') = Tr 



■-u'D^u' — u'ii + uu'Au — u'V {uu) 



(6) 



where V is the differential operator Vijs = i {kj6is + kgSij) /2 (in the Fourier space). The field 
theory with the action functional (jH} has been discussed in details in [H]. The functional 
averages computed with respect to the statistical weight exp 5* can be expanded into the series 
of Feynman diagrams drawn with the interaction vertex in which the spatial derivative V{k) is 




encircled and the arrow denotes the u'-tail, and two propagators (lines) which have the following 
analytical representations (in the Fourier space) 

A - A - ^^^^ (7\ 

The inverse Fourier transform of I\uu' shows that it is retarded, A„u' oc Q(t — to)- Graphs for 
the Green functions in the theory with the action functional © exactly match the well-known 
diagram expansion of Wyld [U] developed by the direct iterations of stochastically driven NS 
with respect to its nonlinear term with the consequent averaging over the Gaussian distributed 
random force. 

One can hardly formulate a general method converting a problem of stochastic dynamics 
into a field theory (there are many, ) but once it is done, it allows for applying the powerful 
techniques within the well controlled approximations. 



3 Branching representations to the Green function of Navier-Stokes 

Many dynamical systems are driven by the non-random external forces. It is worth to 
mention that if one assumes |^| — > (and consequently 0) in the action functional 

it is retained, but then the diagram series is trivial since the velocity field propagator vanishes, 
= 0. However, the diagram series would be recovered by means of regular external forcing. 
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ij + di iujUi) = uAuj + S {x — xq) S (t — to) 



u 



(8) 



supplied with the localized integrable initial condition uq = u (0, x) that means the external 
force / (x, t) = 6 {x — Xq) 6 (t — to) acting on the system. The relevant action functional 



S {u, u') = Tt[—u'u — gvu'V (uu) + uu'Au] + u' (xg, to) 



(9) 



includes the ultra-local interaction term u' (xo,to)- To obtain a formal expansion parameter in 
the perturbation theory for 0, we have inserted the coupling constant gv = 1 m front of the 
interaction term. Despite its physical dimension [g] = —[ly], its conventional dimension would 
be different. The appearance of ultra-local terms in a field theory had been discussed in jHj. 
The basic symmetry of (10) is the Galilean invariance. 



Ua (x, t)-^u{x + Xa (t) , t) - a (t) , u' (x, t) ^ u' {x + {t) , t) 



(10) 



where a [t) is an arbitrary function of time decaying rapidly as \t\ oo, and Xa (t) = J a {T)dT. 



The transformations (fTUI) define the set of orbits in the configuration space along which the 
functional averages do not change. It follows that the functional integral itself is proportional 
to the volume of such orbits. The Green function G{x, t; xq, to) for the Cauchy problem (jH)) can 
be computed as the functional average, 

ff u{x, t) exp S{u, u')DuDu' 



G{x,t;xo,to) = (u) 



-, Sq = Tr[—u'u + z/u'Ati] 
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// exp (So) DuDu' 

in which S is the action functional The above functional average is interpreted as an 

infinite diagram series, 

G(x,t;xo,to) 



(12) 




where the arrows on the bold lines indicate the time direction; the thin arrows display the 
momentum flux traversing the correspondent branch of a graph. The moments are to be 
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conserved at each node of any graph. In comparison to the standard Wyld diagrams the 
graphs sketched in (fT^ contain the local interaction vertex with any number of u'-tails 

that recovers the diagrams even if the only line Au^' is retarded. The similar diagrams 
have been derived for the problem of nonlinear diffusion [Oj, however, in contrast to them, each 
vertex in the diagram series (fT^ contains the differential operator V. The coupling constant 
g = l/u plays the role of an expansion parameter, so that any graph has the amplitude factor 
g''- where i is the number of loops presented in it. 

The diagram expansion for the Green function would have a definite physical meaning if it 
converges. The standard analysis of ultraviolet divergences of graphs is based on the counting 
of relevant canonical dimensions. Dynamical models have two scales, the time scale T and 
the length scale L, consequently the physical dimension of any quantity F can be defined as 
[F] = L~'^fT~^p , in which dp and are the momentum and frequency dimensions of F. 
In diffusion models, these dimensions are always related to each other since dt ~ in the 
diffusion equation that allows us to introduce a combined canonical dimension, dp = dp + 2dp. 
One can check out that each term in Q is dimensionless if the following relations hold: du' = 
0,du = d,di, = —2 + 2 ■ 1 = 0, and dg = 2 — {d + 1) . The field theory (jH)) is logarithmic (the 
conventional dimension of the coupling constant dg = 0) for the Burgers equation {d = 1), while 
in two dimensions dg = —1, and dg = —2 for = 3 (the NS equation). Thus, in the infrared 
region (small moments, large scales) the diagram series in g define just the corrections to the 
diffusion kernel as > 2. However, in the case of Burgers equation, all diagrams look equally 
essential in large scales. 

The diagrams diverge in the ultraviolet region (large moments, small scales) if their canonical 
dimension 

dr = -duNu - du'Nu' > 

where and Nu' are the numbers of corresponding external legs in the graph F. For the 
Green function (lllj). we have Nu' = and Nu = 1. Therefore, there is no any ultraviolet 
divergent graph in the diagram series ()12j). For the first glance, it seems that any graph 
having no external m— legs {N^ = 0) and any number of auxiliary fields u' {N^' > 0) should 
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diverge since du' = 0. However, such a graph is also convergent in small scales because of 
the derivatives V{k) which are always taken outside the graph onto the external m'— legs that 
effectively reduces its canonical dimension to c/p = dr — N^/ < 0. Therefore, the theory Q has 
no ultraviolet divergences and need not to be renormalized. 

Each Feynman diagram in ()12|1 series corresponds to a certain gyro-dynamical process con- 
tributing into G(x, t; Xq, to)- For the convenience of readers, we give the parallel interpretation 
of diagrams in terms of multiplicative branching stochastic processes discussed in |21 ITU] . 
The very first diagram is presents the solution of diffusion equation (in d-dimensional space), 

A(.-.„..-.„)^ -'-(-';-f/^;';-'°'> . (i3) 

{Anu {t -to)) 

It describes the simple viscous dissipation of a vortex with no bifurcations. The second graph 
in (jl2p corresponds to a process of twofold splitting. Under the spatial Fourier transformation, 
it is equivalent to the following analytic expression: 

oo 

Fi (k, t) = -9j dt'V (k) A (k, t - t') j -^A (q, f) no (q) A (k - q, t') (k - q) , (14) 



where A(k, t) is the spatial Fourier transform of the diffusion kernel (14) and mq (k) is the 
Fourier spectrum of initial condition. It is worth to mention that the diagram expansion for 
the Green function can be also discussed for the function defined a the finite domain supplied 
with a periodic boundary condition. In the latter case, it has a discrete set of harmonics, and 
the integral (fT^ turns into sums. 

The time integration in ()14j) can be performed easily, 

rfq r ^^^^t o,,/^_i, .2^ 



ri = -^V(k)exp(-z/A;^(t-to)) / 

J a 



' exp(-2z/(q-k-g ) to) 

uo (q) Mo (k - q) - 



2z/ (q ■ k — g 



2^ 



(15) 

The singularities in Fi appear at q = and q = k, when the vortex does not bifurcate. The 
remaining momentum integral in (fT3|) can be interpreted as an expectation value, 

oo 

Pi (k) = y dry [uo (q) Mo (k - q) exp {-2u (q ■ k - g^) r)] , (16) 



to 



over the Poisson process of vortex bifurcation at momentum k. Consider the statistical ensemble 
of the pitchfork bifurcations (see Fig. 2. a) relevant to the set of possible branching random walks 
discussed in jl]. 

For each simple tree, the processes start at the vortex with the fixed momentum k > 1. The 
branching random walk holds in state k for an exponential time characterized by the parameter 
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uk"^. When clock rings, the vortex bifurcates into two offspring with moments j G [1, A; — 1] and 
k — j. They exist for the exponential times accordingly to their moments and then dissipate 
(that is figured by the crosses in Fig. 2. a). In accordance to the rules of pLj, at each vertex of 
momentum q which dies, one attaches the initial condition UQ{q). The solutions of the fluid 
flows equations (defined in a finite domain with periodic Li initial data) are the expectations of 
a multiplicative functional that is constituted by the recursive application of the following rule 
to each node: —i times the product of the pair of values attached to the corresponding pair of 
input nodes. Each pitchfork bifurcation tree shown in Fig. 2. a presents one of possible primitive 
realizations of the recursive multiplicative algorithm from pp. The analytical expression ()14|) 
for the second diagram sums up the contributions coming from all pitchfork trees in the forest 
shown in Fig. 2. 

The two-loop diagram in (jl2j) represents the forest of multiplicative branching random walks 
with the double pitchfork bifurcations schematically shown on Fig.2.b. The graph corresponds 
to the following analytical expression: 

/\/(k)exp(-z/P(t-to)) xp2(k) (17) 

where the functional 

oo 

P2 (k) = [df [ -^wo (k - p) (p) pi (p) exp (-2Z/ (p ■ k - p') t') 



to 



and pi (k) is defined in (17). The diagrams can be computed provided an integrable initial 
condition 

«ow=Ep;^. (18) 

m 

is given. In ()18|) . the sum is taken over m odd if the space dimension d is even, and over m 
even otherwise. For each term in (jl8|) . one can compute the functional p£ {k) of arbitrary order 



'p ( d m+1 \ 




(19) 



r(^) 

where = id — {i + 1) (m + 1). Accounting for several terms in (fTH|) makes the analytical 
computations very difficult, nevertheless they can be performed numerically. 



4 One-dimensional diffusion with bifurcations 



In the present section we study the correction to the ID diffusion spectrum arisen due to the 
Poisson stochastic process of momentum bifurcations. For the sake of simplicity, we consider 
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the one dimensional diffusion problem supplied with the initial condition in the form of an 
infinite power series, 

uo (k) = y A;-™, A; > (20) 

that corresponds to the function uq (x) = —sin (x)/2 in the real space. The fc-spectrum of the 
unperturbed diffusion kernel at t = 5, z/ = 0.005 is shown by the solid line on Fig. 3. The dash 
line presents the G (/c)-spectrum accounting for the corrections due to the pitchfork bifurcations 
with 

{k) = ^ (cos {k) (Ci {2k) + Ci {-2k) - log (-A;^) - 2 log 2 - 27) + 2Si {2k) sin {k)) (21) 

where Ci {k) is the cosine integral, Si {k) is the sine integral, 7 is the Euler's constant. The 
G {k)- spectrum accounting for the corrections due to both the second and the third diagrams 
in (jl2p is displayed on Fig. 3 by the dash-dot line. We are not assured that the expansion 
parameter g is small, and the expansion (fT^ is the asymptotic series. The fact that the latter 
spectrums coincide gives us only the indirect evidence in favor of relatively fast convergence 
rate of diagram contributions in the series. We discuss this problem in more details in the 
forthcoming section. 

Bifurcation of vortexes is the Poisson stochastic process developing with time. On the 
outlined graphs, we have displayed the exponential decay of G-spectrums at given momentum 
k. In the real space, the Green function G comprises of several typical patterns of sizes oc 1/km 
where km are the local maximums of the spectrum presented on Fig. 3. 

5 The large order asymptotic analysis for the Green function: in- 
stanton approach and Borel summation 

We have mentioned in the previous section that the expansion parameter g in (fT^ is not 
small. To get the information on the convergence of asymptotic series, we study the asymptotic 
behavior of large order coefficients G*-^^ in the diagram series 

N 

g(,) = 5:g(^V, 

k=l 

as iV — > 00, by the asymptotic calculation of the Cauchy integral, 

G(^) = -LI^G (x, t; xo, to) exp {-N log g) , (22) 
2m J g 

in which G{x,t; XQ,to) is the functional integral The contour of integration in ()22|) em- 

braces the point g = in the complex plain. We compute the functional integral by the steepest 
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descent method supposing that is large (instanton method). Instanton approach has been 
apphed to the various problems of stochastic dynamics [HI [T21 HSl UH Hi] • 

Following the traditional instanton analysis, we perform the rescaling of variables in the 
action functional ^ in order to extract their dependence of N, 

u 



uVN, u' uWN, g g/N,u ^ Nu,x ^ xVN. (23) 

This dilatation keeps the action functional Q unchanged, thus each term acquires the multiplier 
N and then formally gets the same order, as the logg in ()22j). The corresponding Jacobians 
from the numerator and the denominator of (jllj) cancel. The saddle point equations are 

ii' + gvuV {u) + vAu = (24) 

u + gvV (uu) - uAu - 6 {x - xq) 6 {t - to) = (25) 

uu'V [uu) = - (26) 
9 

The first equation (j^^ recovers the original Cauchy problem. The additional equations occur 
within the framework of our approach since we analyze the dynamics as a kind of Brownian 
motion. Eq. for the auxiliary field is characterized by the negative viscosity, and therefore 
u'(t > to) = 0. The last equation determines the saddle-point value of g that allows to exclude 
the interaction from the saddle-point equations, 

m'Am = -l+M'5(x-xo)(5(t-to), ^27) 
uA*u' = 1, 

in which we have introduced the diffusion kernel A = dt + I'k'^, and its complex conjugate (in 
the Fourier space) A* = dt — uk"^ . It is interesting to note that the possible solutions of (jTTj) 
should meet the anti-commutation relation, 

u' Au + uA*u' = u'{xo, to)S{'x — xo)5(t — to)- 

Bifurcations of vortexes arisen due to the nonlinearity of hydrodynamic equations do not con- 
clude into a critical regime in the model 0, and therefore the time spectrum in the nonlinear 
model is the same as for the free diffusion equations, r~L^, that is the reason for the branch- 
ing process to be Poisson distributed with the characteristic time l/uk"^. One can see that the 
saddle-point configurations which enjoy ()27|) should not depend upon the Poisson branching, 
therefore, the solutions could exist before the bifurcations commence that is as t — > to- Taking 
into account that 

6 (x) = lim — --, 

e^OTT {x^ + e^) 
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one can find that in the limit t ^ to (j?fjl is satisfied by the following radially symmetric 
solutions, 

u (r, t) = J{r-rof + {t-to)\ u' (r, t) = i^hzAu (r, t) , (28) 
* t — to 

^r., ^ ^a.c.a. (-^) ..'(M, ^ -f <^-'"-;°''"'— (29) 

TT \2z/(t — to)/ |r — ro| 

TT V |r — FqI / 2z/ |r — TqI 



The solutions ()28ll30p are displayed on Figs. 4-5. 

The auxiliary field u' in (j28|l has a pole at t = to, and then from (j26|) it follows that (7* = 
and therefore lays inside the integration contour in (j22I)- In contrast to it, in (jHUj) . u as 
t ^ to and consequently, ^ 00 that is definitely outside the integration contour. Eq. P^jl is 
a subtle point since u{r,t = to) = 9{0)/2, and the position of saddle-point configuration charge 

depends upon the conventional value 6^(0). While estimating the functional Jacobian in 
we had assumed that 6{0) = 0. Following such a convention, one can conclude that in this 
case it is also (7* ^ 00, so that being interested in the large order asymptotic behavior of G^^'' 
we do not need to take the solution ()29|) into account. Even if one takes ^(0) as finite, then 
< g^ < 00, and it is always possible to distort the integration contour in such a way to avoid 
(7* to be encircled. Therefore, is the only solution we need, and below we use 

5f* = hm 5f ~ — -. 

Fields M, u\ and the parameter g fiuctuate around their saddle-point values u*, u'^, and g^,. By 
means of the standard shift of variables, u = + 6u, 6u' = u'^ + 5u' , g = g* + Sg, one makes 
them fiuctuate around zero, so that 6u{oo) = 0, 6u'{oo) = 0. Moreover, Tr 6u = Tr 6u' = 
due to the isotropy of fiuid, and therefore these fiuctuations are posses the central symmetry, 
6u = 6u{r,t), 6u' = 6u'{r,t), the same as the saddle point configuration. 

The contour of integration on the variable 6g now passes through the origin point, and is 
directed there oppositely to the imaginary axis. The integral on 6g is conducted now on a 
rectilinear contour in complex plane {too, —ioo) (in accordance to the standard transformation 
of a contour in the method of the steepest descent). At the turn of the integration contour 
5g i5g, the multiplier (— z) appears so the result G^^^ is real. The contribution to the 
Cauchy integral comes from the pole 6g = —g^ which tends to zero as t ^ to- The values of 
functional integrals on the saddle point configurations are proportional to the entire volume of 
functional integration, they are cancelled in the numerator and denominator simultaneously. 
While calculating the fiuctuation integral, we take into account that the first order contributions 
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in 6u and 6u' are absent because of the saddle-point condition. We also neglect the high-order 
interactions between fluctuations, 0(50^), 0(50^), etc. to arrive at the Gaussian functional 
integral, 



NatN +1/2 f- r T\T 

G^''\r,5t) ~5i^o ^^^^ / / VSu VSu' exp--Tr 



Su'A Su + 6u { -dr 1 Su 



r 

r 



(31) 



Performing the usual rescaling of fluctuation fields 

we compute the Gaussian integral, with respect to 6u first, and then the resulting Gaussian 
integral over 6u', 



G(^)(r,5t) c^st^o iV^-'/'exp (-iVlog^?.) det A'^ (^i + O (^^^^ . 



(32) 



The kernel of operator A^^ is the Green function of the linear diffusion equation (fT^ . Using 
the Stirling's formula, one can check that the coefficients G^^^ of the asymptotic series (|T^ 
demonstrate the factorial grows (like in the most of quantum field theory models ): 



G(^) (r, 6t) ^5*^0 ^ ^ (1 - log 9*) det A"! (^1 + O (^^^ ^ . 



(33) 



Therefore, the asymptotic series (|12j) can be summed by means of Borel's procedure. It consists 
of the following transformation of series ()12|) 

dTg^''\gTfe-^ (34) 

iV iV iV 

where Q^^^ = G^^yT(N + 1) are the new expansion coefficients which do not exhibit the 
factorial growth. We change the orders of summation and integration in ()34|) that is indeed 
incorrect from the mathematical point of view, but is broadly used in physics, 

J2 dTg^''\gTfe-^= dTe--J2^^''^^9rf. (35) 
^ Jo Jo ^ 

We sum over in the rhs of ()35p. 

detA-i r^iar)^ N(i lo^a) detA-^ ^, / g\ 

/ dr e^^ > ^ ^ gjV(i-iogg*) ^ / e^^log 1 - r— , (36) 



and integrate it over r. 



G{g) ^st^o det A-i + i-Ei exp (^-^^^ (37) 
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where Ei(x) is the exponential integral. The asymptotic diffusion kernels accounting for the 
corrections due to the bifurcation processes, 



G{x - xo,t- to) ^t- 



t-tn 



{x - XqY 



exp 



t-to 



{X - XqY 



(38) 



are sketched on Fig. 6. by the dash-dot hues. The unperturbed diffusion kernels are shown by 
the solid hues. 



6 Discussion and Conclusion 

In the present paper, we have studied the hydrodynamical flows of incompressible fluids as 
the Brownian motions over the space of fluid velocity configurations. The essential point of our 
approach is that we have used the field theory formulation of the dynamics which allowed us for 
the implementation of various powerful technics borrowed from the quantum field theory. In 
particular, we have investigated the Green function for the nonlinear hydrodynamical equations. 

In the framework of proposed approach, the Green functions can be calculated in the form of 
well defined asymptotic series. The numerical simulations have demonstrated that these series 
converge relatively fast. We have studied the high order coefficients of the asymptotic series 
around the unperturbed diffusion kernel describing the pure relaxation dynamics. We have 
shown that similarly to the most of quantum field theory models the high order contributions 
into the Green function exhibit a factorial growth. We have summed over the asymptotic series 
by means of the standard Borel summation procedure and found the closed analytical form for 
the corrections (up to an infinite order) to the pure diffusion kernel. 

The value of our paper is that it provides a possible ground for the optimization of existing 
numerical simulation algorithms for the large scale simulations in hydrodynamics. In particular, 
it can be used in studies of the cross-field turbulent transport in plasmas. Moreover, we have 
developed the method that can be successfully used in the turbulence studies. It is worth to 
mention that in the present paper we have constructed the asymptotic solution for the Navier- 
Stocks equation in the "vicinity" of standard diffusion kernel. However, other asymptotical 
regimes would also be possible and have to be studied. 
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Figure 1: The expectations of the multipUcative binary branching process for the 3D FNS forced in the small 
moments k ~ 0{l) (in the large scales). The binary trees have been sampled independently for each velocity 
component, but the total number of samples is limited just to 250 for each point {k,t). The set of randomly 
generated 250 samples obviously does not provide a reliable base for the turbulent statistics for the large k that 
is evident from the excessively strong trembling of data worsening as the momentum k grows up. 
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Figure 2: a). The statistical ensemble of the pitchfork bifurcations relevant to the set of possible branching 
random walks. The branching random walk holds in state k for an exponential time characterized by the 
parameter uk^. When clock rings, the vortex bifurcate into two offspring with moments j G [1, fc — 1] and k — j. 
They exist for the exponential times accordingly to their moments and then dissipate that is figured by the 
crosses, b). The branching random walks with the double pitchfork bifurcations. 




2 4 66 m 1Z 14 16 18 20 



Figure 3: The fc-spectrum of the unperturbed ID diffusion kernel at t = 5, = 0.005 is shown by the solid 
line. The dash line presents the G (fc)-spectrum accounting for the corrections due to the pitchfork bifurcations. 
The G (k)- spectrum accounting for the corrections due to both the second and the third diagrams in H12|l is 
displayed on Fig. 3 by the dash-dot line. On the outlined graphs, we have displayed the exponential decay of 
G-spectrums at given momentum k. 
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Figure 4: The radial profiles (for 5 consequent time steps) of the radially symmetric velocity field u{r, t) and 
the auxiliary field u'{r,t) which meet the saddle-point equations with respect to the high order asymptotic 
contributions into the Green function. These solutions correspond to the effective value 5* — > of the coupling 
constant in the hydro dynamical equation and contribute the asymptotic close to the standard diffusion kernel. 
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Figure 5: The radial profiles (for 5 consequent time steps) of the radially symmetric velocity field u{r,t) and 
the auxiliary field u'{r,t) which meet the saddle-point equations with respect to the high order asymptotic 
contributions into the Green function. The solutions correspond to the effective value gf* — > oo of the coupling 
constant in the hydrodynamical equation. They do not contribute into the high order asymptotic coefficients. 
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Figure 6: The profiles of standard diffusion kernel (the solid line). On the left: G{r) at several consequent time 
steps; on the right: G{t) at several distant points. The dash-dot lines present the asymptotic kernel (as t — > to) 
accounting for all bifurcations (up to an infinite order) of vortexes. 



